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Abstract 

In 1934 J. Leray proposed a regularization of the Navier-Stokes equations 
whose limits were weak solutions of the NSE. Recently, a modification of the 
Leray model, called the Leray-alpha model, has atracted study for turbulent 
flow simulation. One common drawback of Leray type regularizations is their 
low accuracy. Increasing the accuracy of a simulation based on a Leray regu- 
larization requires cutting the averaging radius, i.e., remeshing and resolving 
on finer meshes. This report analyzes a family of Leray type models of ar- 
bitrarily high orders of accuracy for fixed averaging radius. We establish the 
basic theory of the entire family including limiting behavior as the averaging 
radius decreases to zero, (a simple extension of results known for the Leray 
model). We also give a more technically interesting result on the limit as 
the order of the models increases with fixed averaging radius. Because of 
this property, increasing accuracy of the model is potentially cheaper than 
decreasing the averaging radius (or meshwidth) and high order models are 
doubly interesting. 
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1 Introduction 

In 1934 J. Leray |Leray34a| , |Leray34b] studied an interesting regularization of the 
Navier-Stokes equations (NSE). He proved that the regularized NSE has a unique, 
smooth, strong solution and that as the regularization length-scale 5 — > 0, the 
regularized system's solution converges (modulo a subsequence) to a weak solution 
of the Navier-Stokes equations. This model has recently been attracting new interest 
as continuum model upon which large eddy simulation can be based (see, e.g., the 
work of Geurts and Holm |GH03j and Titi and co-workers |CHOT05j . |CTV05j . 
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|VTC05] ). If w denotes a local spacial average of the velocity w associated with 
filter length-scale 5 , the classical Leray model is given by 



(1.1) 



d t w + w • Vw — i/Aw + Vg = f , 
V • w = 0. 



Leray chose w = gs*w, where g$ is a 




averaging radius 5. In (11.11) . v 



is the kinematic viscosity, (w, q) denotes the model's velocity and pressure, fl is the 
flow domain and f is the body force, assumed herein to be smooth and divergence 
free. We take fl = (0, 2tt) 3 , the initial condition 



and impose periodic boundary conditions (with zero mean) on the solution (and all 
problem data) 



The Leray model is easy to solve using standard numerical methods for the 
Navier-Stokes equations, e.g., [LMNR06J, and, properly interpreted, requires no ex- 
tra or ad hoc boundary conditions in the non-periodic case. However, as an LES 
model it has three main shortcomings: 

• The Leray model's solution can suffer an accumulation of energy around the 
cutoff frequency (Geurts and Holm iGH03j). 

• The Gaussian filter is expensive to compute. 

• The model f 1 1.1 1) has low accuracy on the smooth flow components, e.g., | |ujvs£:— 
w|| = 0(5 2 ) , Section 4. 

The accumulation of energy around the cutoff length-scale can be ameliorated 
by new ideas in eddy viscosity which focus its effects on the smallest resolved scales 
or by time relaxation models with similar motivations, |Guerj . |SAK01a] . |SAK01b"] . 
[S AK02J . The expense of computing the filtered velocity is reduced (as proposed by 
Geurts and Holm |GH03j ) if the Gaussian filter is replaced by a differential filter (in- 
troduced into LES by Germano [Ger86j ). The resulting combination of Leray model 
plus differential filter, called the Leray-alpha model, has attracted an explosion of 
recent interest because of its theoretical clarity and computational convenience. For 
recent work on the Leray-alpha model, see Geurts and Holm [GH03j . [GH05] , Guer- 
mond and Prudhomme [GP05] . Cheskidov, Holm, Olson and Titi [CHOT05J, Ilyin, 
Lunasin and Titi [ILT05j . Chepyzhov, Titi and Vishik [OTV05] . |VTG05j (among 
other works). Finally, there is the issue of the accuracy of (II. ip upon the large 
scales / smooth flow components which must be improved if (ll.ip is to evolve from 
a descriptive regularization to a predictive model. 



w(0,x) = uq(x) in fl , 




■"■By other choices of convolution kernel, differential filters, sharp spectral cutoff and top-hat 
filters can be recovered. 



In this paper we complement the above work on (11. II) by presenting the analysis 
of a new (and related) family of Leray-deconvolution models which have arbitrarily 
high orders of accuracy and include the Leray-alpha model as the zeroth order 
(N = 0) case. To present the Leray-deconvolution models, some preliminaries are 
needed. Although any reasonable filter can be used, for definiteness, we fix the 
averaging to be by differential filter. Thus, given 27r-periodic free divergence field w 
with zero mean value, its spacial average over 0(5) length-scales, denoted w, is the 
unique 27r-periodic solution of the Stokes problem^ 

(1.2) - 5 2 Aw + w + Vtt = w inR 3 , V ■ w = 0, / w = 0. 



n 



It can be shown that it is a constant in the equation above and therefore the pressure 
term disappaers. Given a filter, the (ill-posed) deconvolution problem is then: 

given w, find a (useful and stable) approximation to w. 

Let Dn denote the map: w — > chosen approximation to w, introduced in section 
2. The approximate de-convolution operators, iV = 0,l,2,---,we consider have the 
asymptotic accuracy and stability properties (see Section 2) 

(1.3) = DnJ + 0(5 2N+2 ) , for smooth 0, 

(1-4) ||£>JvlkH -»H ) < N+l< OO, 



where 
(1.5) 

The Leray-deconvolution model is then V ■ w = and 



H = jv G Lf 0C (JR 3 ), 2tt - periodic, V • v = 0, j v = oj . 



(1.6) <9tw + (D N w) ■ Vw - t/Aw + Vg = f in R 3 x (0, T), / w = 0, 

Jn 

subject to initial and 2n— periodic boundary conditions. Because of (jl.3p . the formal 
accuracy of (1.4) on the smooth flow components as an approximation of the Navier- 
Stokes equations is 0(5 2N+2 ), N = 0, 1, 2, • • •, Section 4. When N = 0, D w = w so 
(jl.6p reduces to the Leray-alpha regularization. 

The behavior of the model as iV increases for 5 fixed is beyond known Leray- 
type theories and relevant for practical computation. Indeed, decreasing 5 for fixed 
iV requires reducing the computational meshwidth (a process which increases the 
storage and computing time dramatically) and resolving (1.1) or (1.2). On the other 
hand, increasing N for 5 fixed requires only the solution of one additional Poisson 
or Stokes-type problem per deconvolution step. The main theoretical contribution 
of this paper is, in Section [31 to resolve this limiting behavior of the model. We 



2 This filter is very close to the Gaussian filter; it is the first sub-diagonal Pade approximation 
thereof (as in the rational model |GL00j ): in the Camassa-Holm/alpha model [FHTOlj . (|1.2p is 
often called a Helmholz filter. One possible extension of this differential filter to non-periodic 
boundary condition is given in Manica and Merdan MM06J. 



first prove existence and uniqueness of a smooth solution to (11.61) . then we show 
that (modulo a subsequence), for fixed 5, as iV — > oo the model solution w = w(N) 
converges to a weak solution of the Navier-Stokes equations. To our knowledge, this 
is the first result on the limiting behavior of a family of turbulence models as the 
order of accuracy of a family of models on the large scales increases. The difference 
between the two limiting cases is, loosely speaking, that as 5 — > , w — ► w in every 
reasonable sense and the deconvolution process inherits this: Djyw — > w for N 
fixed and as 5 —>■ 0. However, since deconvolution is ill-posed and the deconvolution 
operators Dn are only an asymptotic (in 5 for very smooth solutions and for N 
fixed) inverse, the limiting behavior as iV — > oo is more delicate. 

Another new features include an 0(5 2N+2 ) error bound for the energy norm of the 
model's error (section 0J: ||w(5) — UjvsjsII, provided the underlying solution of the 
Navier-Stokes equations is sufficiently smooth, and an estimate of the time averaged 
error < ||V(uat,s£ — w)||jy >2 both for general weak solutions of the Navier-Stokes 

5 

equations and for weak solutions having the k~s energy spectrum typical observed 
in fully developed turbulent flows. (These first of the two estimates is connected 
to related results for approximate deconvolution models in [LL06a| . [LL06bJ, [LL03J 
and (DE06) and the second is an extension of the authors work in |LL06b| .) 

Of course, ultimately, practical computations require analytic guidance in bal- 
ancing the computational meshwidth with 5 , N and other model and algorithmic 
parameters. 



1.1 Notation and preliminaries 

We use || • || to denote the L 2 (Q) norm and associated operator norm. We always 
impose the zero mean condition J n <fidx = on = w,p, f and w . Recall that 



(1.7) H = jw G Lf 0C (Hr), 2vr - periodic, V ■ w = 0, J w = 
We also define the space function 

(1.8) Hx = |w G H^ 0C {JR 3 ), 2tt - periodic, V • w = 0, J w = oj . 

We shall define the space H s for every s in the same way. We can thus expand the 
velocity in a Fourier series 

w(x, t) = w(k, t)e _lk ' x , where k G Z 3 is the wave-number. 

k 

The Fourier coefficients are given by 

w(k,t) = ^p^w(x,t)e- ik ' x rfx. 

Magnitudes of k are defined by 

\k\ = {\h\ 2 + \k 2 \ 2 + \k 3 \}K 



max{|fci|, \k 2 \, \h\}. 



2tt 

The length-scale of the wave number k is defined by / = — — . Parseval's equality 

k 



oo 



implies that the energy in the flow can be decomposed by wave number as follows. 
For w G H , 

pi p /jiw ( x, t) i^=x:i|w(k, f )p= 

= E( E 5|w(k,«)| 2 ] .where k£Z 3 . 

k \k-K\k\<k J 



Moreover, for s G IR, 

;i.9) 



H s = jv G #f 0C (IR 3 ), 2vr - periodic, V • v = 0, w(0) = 0, |k| 2s |w(k, t)\ 2 < oo 

We define the H s norms by 

(1.10) ||w||^ = ^|k| 2s |w(k,t)| 2 , 



k 

2 



where of course ||w||q = ||w|| 2 . It can be shown that when s is an integer, ||w 
||V s w|| 2 (see |DG95j ). 

1.2 About the filter 

Let w G H and w G Hi be the unique solution to the Stokes problem. 
(1.11) - 5 2 Aw + w + Vtt = w inJR 3 , V ■ w = 0, / w = 0. 



-ik.x 



Jn 

It is usual in deconvolution studies to denote the filtering operation by G so that 
w = Gw. Writing 

k 

it is easily seen that Vn = and 

w(M) 



;i.i2) w(x,t) = J2i 



-ik.x 



- , 5 2 |k| 

k 1 



Then writing w = G(w), we see that in the corresponding spaces of the type H s 
the transfer function of G, denoted by G is the function 

1 



G(k) 



l + 5 2 |k| 2 ' 

and we also can write on H s spaces type 

(1.13) -5 2 Aw + w = w inM 3 , V ■ w = 0, / w = 0. 



n 



Moreover, one notes that the transfer function depends only on the modulus of the 
wave vector k. Therefore, by noting k — |k|, we shall write in the following G(k) 
instead of G(k). 



2 Approximate de-convolution operators 



The de-convolution problem is central in image processing, [BB98j . The basic prob- 
lem in approximate de-convolution is: given w solve approximately for w: 

(2.1) Gw = w. 

Exact de-convolution is typically ill-posed. We consider the van Cittert [BB98J 
approximate de-convolution algorithm and associated operators, introduced into 
LES modeling by Adams, Kleiser and Stolz, e.g., [A"S0l] . [A"S02] . [SA99] . jSAKOlaj . 
|SAK01b] . {SAK02J. For each N = 0,1, ■ ■ ■, = D^w is computed as follows. 

Definition 2.1 (van Cittert approximate de-convolution algorithm) Setw = 
w, 

forn = l,2,---,N-l,_ 
perform w n+1 = w n + {w - Gw n } 
Set Wat = Djvw 

By eliminating the intermediate steps, we find 

TV 

(2.2) D N w = - G T W - 

n=0 

The approximate deconvolution operator Dn is a bounded operator (Lemma 12. ip 
which approximates the unbounded exact deconvolution operator to high asymp- 
totic accuracy 0(5 2N+2 ) on subspaces of smooth functions (Lemma 12.21) . We begin 
by summarizing from [DE06J, [BIL06] a few important, known properties of the 
approximation deconvolution operator . 

Lemma 2.1 [Stability of approximate de- convolution] Let the averaging be defined 
by (11. ^P . Then is a self-adjoint, positive semi-definite operator on Ho and 

\\D N \\ = N+1. 

Lemma 2.2 [Error in approximate de- convolution] Let G be given by the differential 
filter (II. ^p . For any w e H , 

(2.3) w - D N w = (I - (-5 2 A + l)- 1 )^^ 

= (-l)N+l 5 2N+2 A N+l^_ S 2 A + 1 )-(AT+1) W _ 

2.1 Some calculations 

Since we consider the periodic problem, it is insightful visualize the approximate 
de-convolution operators D^ in terms of the transfer function of the operator D^ ■ 
Since these are functions of 5k , where k = |k|, and not 5 or k, it is appropriate to 
record them re-scaled by 

k <- 5k. 



Since G = (— 5 2 A + 1) 1 , we find, after rescaling, G(k) = r-r 

1 + k z 

transfer function of the first three deconvolution operators are 



. With that, the 



D 



2k 2 + 1 



1 + 



k 2 + 1 k 2 + 1 ' 

7.2 / 1,2 \ 2 

+ 



and 



k 2 + l \k 2 + l 



These three are plotted together with the transfer function of exact de-convolution 
(k 2 + 1) (in bold). 

Figure 1: Exact and approximate (N = 0, 1, 2)) de-convolution operators 



transfer junctions 




.7 -6 



wave number k 



More generally, we find 

n=0 \ ~ / 

The large scales are associated with the wave numbers near zero (i.e., \k\ small). 
Thus, the fact that D N is a very accurate solution of the de-convolution problem 
for the large scales is reflected in the above graph in that the transfer functions 
have high order contact near k = 0. The key observation in studying asymptotics of 
the model as N —>■ oo is that (loosely speaking) the region of high accuracy grows 
(slowly) as well as N increases. 

The regularization in the nonlinearity involves a special combination of averaging 
and deconvolution that we shall denote by H^. 

Definition 2.2 The truncation operator : H — * H is defined by 



1 + k 2 







N+l 



Hjyw := .Datw = (Dn ° G)w. 



Proposition 2.1 For each N — 0, 1, • • • , the operator is positive semi- definite 
on H . The operator H N : H — > H is a compact operator. Moreover, H N maps 
continuously H onto H 2 . Further, 



i^ArW = -DatW, 

\\H N \\c(H ^H ) = 1 , for each N > 0, 
H-ffjvlk(Hi-Hi) = 1 ; for each N > 0. 

Proof These properties are easily read off from the transfer function Hj^(k) of 
if TV which we give next. For example, compactness follows since \H]y(k)\ — > as 
|&| — > oo. 

Remark 2.1 Similarly, by the same proof, maps H s into itself (see Lemma l275\ 



below), and is compact. Since H N (k) 



< 1, one has for each s > 



(2.4) ||#atw|| s < ||w|| s . 

The Fourier coefficients/transfer function of the operator are similarly easily 
calculated to be (after rescaling by k <— 5k) 



(2.5) H N (k) = 1 



k 2 x N+1 



1 + k 2 

They are plotted below for a few values of N. 



Figure 2: Transfer function of H N , N = 0, 10, 50) 




These plots are representative of the behavior of the whole family. Examining the 
above graphs, we observe that H N {u) is very close to u for the low frequencies/largest 



solution scales and that H^{u) attenuates small scales/high frequencies. The break- 
point between the low frequencies and high frequencies is somewhat arbitrary. The 
following (from [LN06aj ) is convenient for our purposes and fits our intuition of an 
approximate spectral cutoff operator. We take for k c the frequency for which 
most closely attains the value |. 

Definition 2.3 (Cutoff- Frequency) The cutoff frequency of is 



k c := greatest integer 



(MS))- 



In other words, the frequency for which Hn most closely attains the value ~. 

From the above explicit formulas, it is easy to verify that the cutoff frequency 
grows to infinity slowly as iV — * oo for fixed 5 and as 5 — > for fixed N. Other prop- 
erties (whose proofs are simple calculations) of the operator Hn(-) follow similarly 



easily from H N (k) 



k 2 x N+1 



l + k\ 

Lemma 2.3 For all N > 0, H^{k) has the following properties: 

< H N (k) < 1, 
\HN(k)\ — > as \k\ — > oo, 
H N (0) = 1. 

Let k c be cutoff frequency, then 

\k c \ — > oo as N — > oo for fixed S and as 5 — > for N fixed. 

For any fixed value of k (or bounded set of values of k ) 

i^jv(fc) — > 1 as N —y oo for 5 fixed and as 5 —>■ for N fixed 

Proof. This follows from (I2.5p . 

To extract strong convergence of deconvolution, some restriction to the large 
scales is needed (as shown clearly in the above two figures). One way to do this is 
to consider the action of deconvolution operators on trigonometric polynomials. 



Definition 2.4 The space IP n denotes the 2n-periodic, divergence- free, wit hzero 
mean value, vector, trigonometric polynomials of degree < n and U n : Ho — > IP„ 
is the associated orthogonal projection. 



The model error is driven by the error in approximate deconvolution. With this 
definition a critical convergence result on approximate deconvolution is possible. 



Proposition 2.2 Hn is a compact, symmetric positive semi-definite operator. For 
n fixed, H N is strictly positive definite on IP n and 

Hn — > I on TP n in the strong operator norm on Hq. 

Given n there is an N (n) large enough such that 

||n n v|| < 2||#Jv|| , for all N > N {n), 

Further, 

i 

||v|| < ||^v|| < 2||v|| for all v £ P n . 

Proof. This follows from (12.51) . For example, Compactness follows since 
H^ik) — > as k — > oo. 

Other properties of the operator H^(-) follow similarly easily from its transfer 
function. 

Proposition 2.3 Let IT/^ denote the orthogonal H projection into 
span{e* kx : |k| < k c }. For all w £ H .' 

(2.6) (H N w,w) >C\\U k w|| 2 , 

(w-H N w,w) > C\\(I-Tl kc )w\\ 2 . 

Proof The claims follow from the definition of the cutoff frequency, the explicit 
formula for the transfer function and direct calculation. 

The following properties are important to the analysis that follows. 

Lemma 2.4 Let s £ 1R ; w £ H s . Then H N (w) £ H s+2 . 

Proof. Since ifjv( w ) satisfies periodic boundary conditions, it is divergence free 
and has zero mean value by construction. We check the regularity. It is easy checked 
by using a Taylor expansion that there exists a constant C = C(5, N) such that 

(2.7) Vk £ Z 3 , k ^ 0, < H N (k) < C ^' N ^ , with k = k. 

k l 

Therefore, for w £ H s , w(x) = Xlkez w(k)e~ lk ' x 

(2.8) (H N (k)) 2 k 2( - s+2 ^ | w(k) | 2 < C (5, N) 2 k 2s \w(k)\ 2 . 

Recall that here w(0) = 0. Therefore in (12.81) . k ^ 0. Hence, by definition 
(2-9) ||^(w)|| Hs+2 <C(5,iV)||w|| Hs , 

and the proof is finished. 

Remark 2.2 The constant C(6,N) in the above proof blows up as 6 — > or as 
N -> oo. 



Lemma 2.5 1) The operator maps continuously H s into H s and 

(2-10) ||#iv|U(H s ,H s ) = 1- 

2) For every w G H s; H N (w) converges strongly to w m H s w/ien iV — > oo. 
ifjv commutes with the gradient operator: VHn = H^V . 

Proof. 1) Let w = (w 1 , w 2 , w 3 ) G H s . One has 

||w|| 2 =]>> 2 ^(k)| 2 . 
kez 3 

Let us note w(k) = (w^, w 2 ,, w 3 .), H N (w) = (/^(w) 1 , H N (w) 2 , H N (w) 3 ). Therefore, 
the j th component of the vector field H N (w) is given by 

(2.11) H N (wY = H N (k)wi^. 

kez 3 

Consequently 

V • H N (w) = iHjvi^kjwie 1 ^. 
kez 3 

When w G Hi, Vk G Z 3 , kjW J k = 0. Then V • H N {w) = 0. Now recall that 

p \ N+l 

H N (k) = 1 - ' " 



1 + k 2 

One deduces that \H N (k)\ < 1 and 

k 2s \H N (k)w(k)\ 2 < k 2s \w(k)\ 2 
yielding ||i?jv(w)|| s < ||w|| s . Hence 

H N (w) G H s and 1 1 H N ( w) 1 1 Hs < ||w|| Hs , 
showing ||i?Ar||£(H s ,H s ) < 1- Now let k = /c 2 ,0) G Z, k\ ^ 0, k = |k| and define 

* (k) = wrmm (h - (hm Wk(x) = Mk)e ' kx - 

By construction w k G H s , ||w k || Hs = 1, and one also has 1 1 i^";v w k I |h s = H N {k) — > 1 
as k — > oo, making sure that ||.Hjv|U(h s ,h s ) — 1- 
2) Let w G H s . One has 

/ L.2 \ 2(JV+1) 

||w-^w|| s =^ — - A;-|w(k)| 2 
kez 3 v ' 

Let e > 0. Since w G H s , there exists K > such that 



5> 2 1w(k)| 2 <^, 



yielding 

||W - H N W\\ S < ( ; — ^ 1 l|w|| s + 



K 2 \ 2{N+1 \, , e 



1 + K 2 J 11 2 

Now, as K > 0, i^ 2 /(l + K 2 ) < 1, there exists iVo (which depends upon w) such 
that ViV > JV , 

K 2 \ 2{7V+1) ,, e 



w L < -. 



1 + if2 y ii n»- 2 - 

Then V N > iV , | |w — iJ^wl |h s < £■ This shows the strong H s convergence of H^w 
to w as iV ^ oo. 

3) The fact that commutes with V is obvious. 

Remark 2.3 All the results in Lemma 2.5 above are uniform in 5. Moreover, it 
easy checked that for a fixed N, the same convergence result holds when 5 goes to 
zero. 

Remark 2.4 Let w G L 2 ([0, T], H s ). It easy checked by using the same proof as in 
Lemma [2~5\ combined with the Lebesgue monotone convergence theorem that H^(w) 
converges to w in L 2 ([0,T],H s ) as N goes to infinity. 

3 The Leray Deconvolution Model 

3.1 Existence result 

The theory of Leray-deconvolution model begins, like the Leray theory of the Navier- 
Stokes equations, with a clear global energy balance, and existence and uniqueness 
of solutions. 

Let u G Ho, / G H_i. For 5 > 0, let the averaging be defined by ( 11. lip . The 
problem we consider is the following, for a fixed T > 0, find (w, q) 

( wGL 2 ([0,T],H 1 )nL oo ([0,T],H ), d t w G L 2 ([0, T], H_x) 

I gGL 2 ([0,T],L 2 ei . ), 

{ 1 | d t w + (H N (w)V)w - v Aw + Vq = H N (f) in V'([0, T] x R, 3 ), 

[ w(x, 0) = H N (u ) = w . 

where L per denotes the scalar fields in L 2 oc (IR 3 ), 2-7r-periodic with zero mean value. 

Theorem 3.1 The problem (\3.1\i admits a unique solution (w, q), where w satisfies 
the energy equality 

(3.2) \\\™{t)f + v [ [ \Vw\ 2 dxdt' = l-\\H N (u \\ 2 + [ [ H N (f).wdxdt'. 
2 Jo Jn 2 J J n 

Moreover, w G L°°([0, T], Hi) n L 2 ([0, T], H 2 ) 



Proof. For simplification, one notes || ■ \\ Pj g the norm in L p ([0,T], B) and when 
B = H s one shall note 1 1 • | | P)S for the simplicity. One also considers the space 



(3.3) V = L 2 ([0,T],H!) nL°°([0,r],Ho), ||w|| v = ||w|| 2jl + ||w| 

Recall that (see [Le06] ). 



oo.Q- 



(3.4) VwGV, Wg[2,6], ||w|| 4r £r < C[|w| 

We introduce the operators: 

r-T 



A(w,v)=v [ [ Vw.Vv, 

B N (w,v) = / / (iJ A r(w)V)w.v = — / / H N (w) <g> w : Vv. 



'o ./o ./n 

We first notice that 

(3.6) |A(w,v)| < f| |w| | 2> i 1 1 v| | 2 ,i- 

Now assume that w G V, v G L 2 ([0, T], Hi). Using lemma [23 combined with an 
integration by parts and the Sobolev Theorem, it is easy seen that there exits a 
constant C which depends on N and such that one has 

(3.7) |5jv(w, v)| < C||w||fr||v||2,i, Bjv(w, w) = 0. 

Here we have remarked that because H 2 C (L°°(IR)) 3 (we are working in a 3D 
configuration), w G L°°([0, T], H ) yields H N (w) G L°°([0, T], H 2 ) (see §ZM)) and 
therefore, 

(3.8) H N (w) GL°°([0,T] xM), H^w)!^ < C7||w|| v 
for a constant C which depends on iV. We now note 

(3.9) W = {w G V, d f w G L 2 ([0, T], H_i)}. 

Notice that any w G W is almost every- where equal to a function in C ([0, T], Ho) 
(see in |RT84j ). Therefore, passing to a quotient space by keeping the same nota- 
tions, we have 

WCC°([0,T],H ). 

Finally when f G L 2 ([0, T], H_i), then H N (f) G L 2 ([0,T],Hi) and the variational 
formulation of Problem (13.11) is: 

Find w G W such that w(0,x)=w , and Vv G L 2 ([0, T], Hi), 
( ' , ' 1U ' ) ^ < <9 t w,v > +B N (w,v) + A(w,v) = / [ H N (f).vdxdt 



o ./n 



Since (13. 6p and (13. 7p are satisfied, the existence of a solution to Problem (I3.10P can 
be derived thanks to the Galerkin method. This is classical and the reader can look 
at |LLe06] and references inside or also in [RT84] for a detailed descrition of such 



kind of proof. By taking w = v as test vector field in ( 13.101) (a legal operation 
here), one directly gets the energy estimate (13.21) . Notice that thanks to the energy 
equality and ||ffjv(uo)|| < H u o||> I l-^iv(f) I l-i < ||f||-i> w does satisfy the following 
estimates: 

(3.11) I|W||L, < l|Wo|| 2 + ^||f||2,-l, ||W||^ < Illwol^ + ^Hflll,.! 

in other words, 

(3-12) ||w|| v <C(i/,||uo|U|f||2,-i). 

Here the bound does not depend on N. Notice also that we can derive an estimate for 
<9 4 w in the space L 2 ([0, T], H_i). However the bound for d t w in this space depends 
on N. 

The pressure q is recovered thanks the De Rham Theorem and its regularity 
results from the fact that Vg G L 2 ([0, T], H_i). 

We now check regularity. Let Dw and Dq be a differential of w and q. Thanks 
to the periodic boundary conditions, one has 

(3.13) d t Dw + (H N (w)V)Dw - uADw + VDq = DH N (f) - (DH N (w)V).w 

One notes that DH N (f) G L 2 ([0, T], (H ) 3 ), DH N (w) G L°°([0, T], (Hi) 3 ). Since 
w G L 2 ([0, T], H ), by the Sobolev imbedding Theorem, (Di7jv(w)V).w is periodic 
and in the space L 2 ([0,T], (L 3 / 2 ) 9 ). Since 1/6 + 2/3 = 5/6 < 1, one has 

(Mjv(w)V).w G L 2 ([0,T], (L 3 / 2 ) 9 ) C L 2 ([0, T], H_x). 

Then the equation (13.131) admits a unique solution in the space V 3 . By using 
the same technique as in |LLe06j . it is easy to show that this solution is equal to 
Dw showing that w e L°°([0, T], Hi) n L 2 ([0, T], H 2 ). In particular, there is some 
constant C (which depends on N) such that 

(3-14) l|w|U,i + | |w|| 2)2 < C. 

It remains to prove the uniqueness. Let (wi, qi) and (w 2 ,g 2 ) be two solutions, 
5w = wi — W2, 8q = qi — q\. Then one has 

(3.15) d t 5w + (H N (w 1 )V)5w - uA5w + \75q = (H N (5w)V)w 2 , 

and 5w = at initial time. All the terms in the equation above being in L 2 ([0, T], H_i), 
one can take Sw G W C L 2 ([0,T],Hi) as test, a legal operation. Since H^(wi) is 
divergence free, one has 



I [ (H N (w 1 )V)5w.5w = 0. 
Jo Jn 



Therefore, 



(3.16) ^ / |()w|- + /' / |V()w| L '= / (H x (,)w)Y)w 2 .<)w 



One has by a part integration, 

J (Hn(5w) V)w 2 .£w = — J Hn(5w) (g) w 2 : V<5w 
By Young inequality, 

J (Hn(5w) V)w2-5w = — A i/Ar(5w)®W 2 : V5w < — / |V<5w| 2 H A |-£/jv((5w)| 2 |w 2 | 

Hence, using (13.141) 

(3.17) A / iJwf + | /" | W w| 2 < ^||w 2 ||^||J^ iV (^w)|| 2 

•J & l Oil 

Therefore, 

^||HI 2 <c(*)||HI 2 , 

where C(t) = — ||w 2 || 2 X) £ L 1 ([0,T]). We conclude that 5w = thanks to Gron- 
wall's Lemma. 

Remark 3.1 It easy checked that the solution (w, g) issuchthatw £ C°([0,T], Hi). 

Remark 3.2 .Efy iterating the process obove, it is easy to prove that when f £ 
C°°([0,T] x R) and u £ C°°(R) f&ot/i space periodic), then w £ C°°([0,T] x R). 

Remark 3.3 The Galerkin approximations to the solution w are under the form 

w n = w„(t,k)e ik - x , w n (0,x) = n n (w )(x). 

|fc|<n 

where the vector (w n (i, — n), w n (i, n)) zs a solution of an ODE and is of class 
C 1 on [0, T] according to the Cauchy Lipchitz Theorem. We know on one hand that 
the sequence (w n ) ngN converges strongly to w in L 2 ([0, T), H ) whether the sequence 
(<9iW„) ng ]N converges weakly to d t w in L 2 ([0, T), H_i) . On the other hand, for any 
smooth periodic field 4> with V • <p = and 4>{T, x) = 0, since w n £ C 1 ([0, T],JP n ), 
a legal part integration yields 

< d t w n , 4>>= / d t w n (t, x). 4>(t, x) dxdt = 



o Jo Jn ^ 

T 



Il n (w )(x). 0(0, x) dx - / / w„(t, x). x) dxcft. 

n Jo in 



Passing to the limit when n — > oo yields 



(3.18) / <9i\v,0>= / w o (x).0(O,x)rfx- / / w(t,x).dt4>(t,x)dxdt. 
Jo Jn Jo Jn 



3.2 Limiting behavior of the Leray-deconvolution model 



The models considered are intended as approximations to the Navier-Stokes equa- 
tions. Thus the limiting behavior of the models solution is of primary interest. There 
are two natural limits: 5 — > and N —>■ oo. The first is the normal analytic question 
of turbulence modeling and considered already by J. Leray in 1934. The question of 
the behavior of the model's solution as iV — > oo is much more unclear, however. 

In practical computations, cutting 5 means re-meshing and increasing the mem- 
ory and run time requirements greatly while increasing N simply means solving one 
more shifted Poisson problem per deconvolution step. Thus, increasing the accuracy 
of the model is much easier than decreasing its resolution. On the other hand, the 
van Cittert deconvolution procedure itself is an asymptotic approximation rather 
than a convergent one and has a large error at smaller length scales. 

Before developing these results some preliminary definitions are needed. 

The notion of weak solution is due to J. Leray |Leray34a| who called them tur- 
bulent solutions. Recall that the Navier-Stokes equations are the following: 

d t u + u ■ Vu — i/Au + Vp = f , 
(3.19) V-u = 0, 

u(0,x) = u 

Definition 3.1 [Weak solutions of Navier-Stokes Equations] Let u G H f G 
L 2 ([0, T], H_i). A measurable vector field u(i, x) : Q x [0, T] — > R 3 is a weak 
solution to the Navier-Stokes equations if 

(i) u G V (where V is defined by ( TO) ), 

(ii) u satisfies the integral relation: 



[u.d t <f> + v[ /vu:V0+/ /(uV)u.0 

ii in Jo Jn Jo T Jn 



<f,<f>> dt)- / u .0(0 



(3.20) 



for all d) G (C°°([0,T] x 1R)) 3 ; space periodic with 0(T,x) = forall x G H 3 and 
such that V • (f> = 0. 

(Hi) [Leray's inequality /the energy inequality] for any t G [0, T] 

(3.21) i||u(f, -)|| 2 + v [ \\Vu(t', -)\\ 2 dt' < h\u \\ 2 + [ < f,u > dt', 

* Jo * Jo 

(iv) lim I \u(t, ■) — unl I = 0. 

t^o 

A N being given, we denote by (w^,^) the unique solution to the Leray- 
Deconvolution problem ( 13. ip . wat(0,x) = w ,jv(x) = ifjv(uci)(x). We prove the 
following 

Theorem 3.2 There exists a sequence (Nj)j^ be such that (wjv-)jeK converges 
to a weak solution u to the Navier-Stokes Equations. The convergence is weak in 
L 2 ([0,T],H!) and strong m L 2 ([0,T],H ) 



Proof. Thanks to the bound (13.111) . from the sequence (wjv) ngn one can extract 
a subsequence (wjv.)jGN which converges to some u, weakly in L 2 ([0, T], Hi). In the 
following, we shall denote by (wj) Jg]N this subsequence and we have to show that u 
is a weak solution to the Navier-Stokes equations as defined in Definition (13. ip . 

Thanks to (13.41) combined with (13. lip , the sequence (Wj)j- e N is a lso bounded 
in L 8 / 3 ([0,T],(L 4 ) 3 ), while (^.(w,-)),^ is bounded in L°°([0, T], H ). Hence the 
sequence (H N .(wj) <8> Wj^-gN is bounded in L 8 / 3 ([0, T], (L 4 / 3 ) 9 ) making B N .(wj,-) 
bounded in the space L 8 / 5 ([0, T], H_ 2 ). 

Obviously, (A(wj, -))ieisr is bounded in L 2 ([0, T], H_x) as well as (H N .(f))j £N 
Writing the equation for Wj under the form, 

(3.22) d t Wj = -A{wj, ■) - B Nj { Wj , ■) + H N .(f), 

one deduces that the sequence (d t Wj)j e -^ is bounded in L 8//5 ([0, T], H_ 2 ). Since one 
has Hi C H C H_ 2 , the first injection being continuous compact and dense, the 
second being continuous and dense, one deduces from Aubin-Lions Lemma (see in 
[JS87]) that the sequence (wj) jeI ^ is compact in L 8 / 3 ([0, T], H ), and by very classical 
arguments in all L p ([0, T], H ), p < oo and also in L q ([0, T], (L 4 ) 3 ) for q < 8/3. 

Even extracting an other subsequence still denoted by (wj) je K, the sequence 
(wj)j G ]N converges to u almost everywhere (use the Lebesgue inverse Theorem). 
Notice that thanks to Fatou's Lemma, for almost every t G [0,T], one has 

(3.23) ||u(t,-)|| < liminf | Iw^t, Oil- 
Then by (ETTTll . u e V. 

We have to show that u satisfies points (i) and (ii) of Definition 13.11 to prove 
that it is a weak solution to the Navier-Stokes equations, knowing already that it 
satisfies point (i). 

We check point (ii). Let (/> G C°°([0,T] x H), space periodic with 0(T,x) = 
for all x G 1R 3 and V ■ 4> = (we denote by T the space made of such fields). Notice 
that (j> G L 2 ([0,T],Hi) and therefore can be used as test vector field in formulation 
(13.101) . It is obvious that one has 

(3.24) lim A{wj, 0) = A(u, 0) = / / Vu(t, x) : V<f>(t, x) dxdt. 

'J^°° Jo Jn 

By using the result in Remark 12.41 it is also obvious that 

(3.25) lim / / H N (f)(t,x) .w j (t,x)dxdt = [ <f,u>, 
3-^°° Jo Jn ' Jo 

as well as thanks to ( 13.181) and since (i?Ar(uo))jvGN converges to u in the space H 
(Lemma 12.51) . 

(3.26) lim / < Wj, 4> >= / u o( x ) • 0(0, x) dx. — / / u(t, x) . d t <f>(t, x) rfxrft. 
' x Jt\ ' Jn ' Jn Jn 



It remains to pass to the limit in the non linearity. Notice that 

\\H Nj {-Wj)-u\\ 2 ,o < \\H Nj (wj- u)\\ 2;0 +\\H Nj (u)-u\\ 2 ,o < 1 1 — u| 1 2,0+ 1 l-^Tzv,- W — "I | a, 

Hence it is easy deduced that (ifjv-( w j))jeiNr converges towards u in L 2 ([0,T],H ) 
and consequently (Hjyjwj) <8> Wj)j £ n converges towards u <g> u in L 1 ([0, T], (L 1 ) 9 ). 
Therefore, 

(3.27) lim B(w,-,<£) = B(u,<f>) = - I [u®u:V<t>=[ /(uV)u.</>. 

" Jo Jn Jo Jn 

Combining (13T24D . fl3T25|) . (13T26D and (l3T2Tj) makes sure that (13T20D is satisfied. 

We now check point (iii). We already know that (9 t w) jGN is bounded in the space 
L 8 / 3 ([0, T], H_ 2 ). Thus, up to a subsequence, it converges weakly in this space to 
some g. Passing to the limit in ( 13. 18ft for <j> G T , one sees that g satisfies 



(3.28) / <g,<f>>= [ u o (x).0(O,x)dx- ! f u(t, x). d t (j)(t, x) dxdt 
Jo Jn Jo Jn 

Hence the following relation holds in the space L 8/,3 ([0, T], H_ 2 ) for all t G [0, T]: 

(3.29) u(t,-) = u (t,-)+ / g{s)ds 

Jo 

Consequently, u G C°([0, T], H_ 2 ). Since H C H_ 2 , the injection being dense, and 
u G L°°([0, T], H ), u(t, ■) does exists for each t G [0,T] and is weakly continuous 
from [0, T] into Ho- Moreover, one has u(0, •) = Uo- Now by weak convergence of 
(wj)j g N to u in L 2 ([0, T], Hi), one has 



(3.30) / / |Vuf<liminf / / |Vw 

'0 Jn 'J^°° Jo Jn 



t p 

2 / / '^™f. 



It is easy checked that (H^^.w^j^ converges to f), u G L 1 ([0,T] x fl), and 
we have previously shown that (H^.(uo))jei<i converges towards u . Since each w 3 - 
satisfies the energy equality ( 13.21) . ( 13.231) combined with ( I3.30P ensures that the 
energy inequality ( I3.21j) holds for almost every t in [0,T]. We have to prove that it 
holds for every t G [0, T}. Let t G [0, T] and (t/Jfcgw be a sequence that converges to 
t and that satisfies (13.211) for each k. We already know that (u(i fc , -))fc e N converges 
weakly to u(t, •) in the space H . Therefore, 

||u(t,-)|| <liminf ||u(t fc ,-)||. 
One deduces from this fact that (13.211) is satisfied because 

[ S \\Vu(s' r )\\ 2 ds', f <f(s , ),u(s',-) >ds'. 

Jo Jo 

are continous functions of s. 



We finish by checking point (iv). We already knows that u(t, •) converges weakly 
to u in H . Therefore 

1 1 Tin 1 1 < liminf \\u(t, Oil. 

t-»o 

As ||Vu(£, -)|| 2 and < f , u > are both in L 1 ([0,T]), one deduces from the energy 
inequality ( 13.211) that 

limsup ||u(£, -)|| < ||uo|| 

t-*o 

Therefore, lim \\u(t, -)|| = ||uq|| which combined to the weak convergence in Ho 
garanties the strong convergence, in particular lim \ \u(t, ■) — u || =0 and the proof 
is finished. 

Remark 3.4 By using same arguments as above, one can show that for a fixed N, 
the sequence of solution to ( Iff, jp converges to a weak solution of the Navier-Stokes 
solution when 5 goes to zero. The proof is left to the reader. 

4 Accuracy of Leray and Leray-deconvolution mod- 
els 

The accuracy of a regularization model as 6 — > is typically studied in two ways. 
The first, called a posteriori analysis in turbulence model validation, is to obtain 
via direct numerical simulation (or from a DNS database) a "truth" solution of the 
Navier-Stokes equations, then, to solve the model numerically for varying values of 

5 and compute directly various modeling errors, such as u — w and u — w . The 
second approach, known as a priori analysis in turbulence model validation studies 
(and is exactly an experimental estimation of a model's consistency error), is to 
compute the residual of the true solution of the Navier-Stokes equations (obtained 
from a DNS database) in the model. For example, to assess the consistency error of 
the Leray (and Leray-alpha model) model, the Navier-Stokes equations is rewritten 
to make the Leray model appear on the LHS as 

a t u + u-Vu- z/Au + Vp - f = V • [uu - uu] in [0, T] x Q. 

The Leray- model's consistency error tensor is then r£ eraH (u, u) := uu — uu . Analy- 
sis of the modeling error in various deconvolution models, various norms and diverse 
settings in |LL06a| . |LL06b| . [BIL06| and |DE06| has shown that the energy norm 
of the model error, ||ujvsb — w modei || or ||ujv5£; — w mod ^||| as appropriate, is driven 
by the consistency error tensor r rather than V - t. 

Thus, an analysis of a model's consistency error analysis evaluates ||uu — uu||. 
In the analysis of consistency errors, there are three interesting and important cases. 
Naturally, the case where u is a general, weak solution of the Navier-Stokes equations 
is most interesting and equally naturally nothing can be expected within current 
mathematical techniques beyond very weak convergence to zero, possibly modulo a 
subsequence. Next is the case of smooth solutions (the classical case for evaluating 
consistency errors analytically). The case of smooth solutions is important for tran- 
sitional flows and regions in non- homogeneous turbulence and it is an important 



analytical check that the LES model is very close to the Navier-Stokes equations 
on the large scales. The third case (introduced in [LL06b| ) is to study time aver- 
aged consistency errors using the intermediate regularity observed in typical time 
averaged turbulent velocities. 

For the Gaussian filter it is known (e.g., Chapter 1 in [BIL06] ) that for smooth 
cf) , <p — (gs * <p) = 0(5 2 ) so that the Leray model's consistency error is second order 
accurate in S on the smooth velocity components: r Leray = (u — u)u = 0(5 2 ). This 
simple calculation shows that the consistency error is dominated by the error in the 
regularization of the convecting velocity. Thus, improving the accuracy of a Leray- 
type regularization model hinges on improving the accuracy of the regularization. 
For the differential filter, from (1.2), <fi — = 5 2 (— A0) so the consistency error of 
the Leray-alpha model (Geurts and Holm [GH03J) is also 0(5 2 ). To make this more 
precise, we begin with a simple lemma of [LL06aj . given in the case of scalar fields 
for the simplicity, the same result being true in the case of smooth vector fields. (We 
include a short proof of this lemma for completeness.) 

Lemma 4.1 Let u = (—5 2 A + l) _1 n. Then for any derivative d^u/dx 13 , with 
multi-index > 



[u — u 



<5 2 ||A— — tu\\ < <5 2 ||A— — rttll , and 



Q X P Q X P Q x 

Hfe?(«-«)II^2 I|V &? U|1 - 

Proof The first follows from the definition of averaging equation, stability of 
averaging and the fact that derivatives commute under periodic boundary conditions. 
For the second, note that the error equation satisfies the equation 

-5 2 A(u -u) + (u-u) = -5 2 Au , 

which is an identity. Differentiating through this equation shows that the derivatives 
of the error also satisfy the same equation (with the derivative of u on the RHS 
instead of u ). Multiplying by (u — u), integrating over f2, integrating by parts and 
using the Cauchy Schwarz inequality gives 

TTMI2 , 11/.. rrr\\\2 / X 2 



^||V(u-u)|r + ||(u-u)||^<^||Vu||||V(u-u)|| 

< 5 2 \\V(u-u)\\ 2 + ^\\Vu\\ 2 , 

and the result follows for the the error (u — u). Since the derivatives of the error sat- 
isfy the same equations as the error, the same proof shows works for the derivatives 
of the error as well. 



4.1 Consistency error of the Leray and Leray-alpha model 



The Leray /Leray-alpha model is the case N = in the family of Leray-deconvolution 
models so we shall denote the consistency error tensor of the Leray-alpha model by 



7"o(u, u) where 

r (u, u) := u (g) u — u g) u , 
where u = (— <5 2 A + l) _1 u (for Leray-alpha), 
and u = gs * u (for Leray). 

For the simplicity, we shall write u Cg> u — u ® u = uu — uu. 

Using estimates of both filter's accuracy, a sharp estimate of the consistency 
error of both can be given. 

Proposition 4.1 The consistency error of the Leray model and Leray-alpha model 
satisfy, for smooth u, 

[ |r (u,u)|dx < C5||u||min{||Vu||,5||Au||}. 
Jq 

Proof By the Cauchy-Schwarz inequality 

|to(u, u)|<ix < | |u — u| 1 1 |u| |. 

n 

Both filters satisfy 

||u-u|| < C5min{||Vu||,5||Au||}, 

e.g., for the Leray-alpha model C = 1, from which the result follows. 

Naturally, for a general weak solution of the Navier-Stokes equations, the indi- 
cated norms on the above RHS may or may not be finite at specific times. Time 
averaged values, however, are always well defined and can be estimated in terms of 
model parameters and the Reynolds number as in [LL06b|. 

Definition 4.1 Let < ■ > denote long time averaging, given by 

1 ' T 



(4.1) < <P >:= LIM^- / <p(t)dt, 

1 Jo 

where LIM denotes the generalized limit for bounded functions introduced in IFMRTOlf 

Definition 4.2 Let e(u) and e denote respectively the energy dissipation rate of 
the flow (the unknown true solution of the Navier-Stokes equations) and its time 
averaged value, defined by 



s(u)(t) — — / u\Vu\ 2 dx , and 
L 6 Jn 

e =< e(u)(t) > . 



Lemma 4.2 Iff 6 L°°([0, T], H_i) then < e(u) >< oo and < e(u) >< oo. 



Proof The NSE result is standard, e.g., Doering and Gibbon [DG95]. The result 
for the model is proven the same way (divide the energy equality by T and take the 
limit as T — > oo. Using Gronwall's inequality, |;||u(T)|| 2 — > 0. The remainder 
follows from the Cauchy-Schwarz inequality). 

Using the estimate of the regularization's accuracy, a sharp estimate of the Leray 
model's consistency error can be given. Naturally, for a general weak solution of the 
Navier-Stokes equations, the indicated norms on the RHS may or may not be finite 
at specific times. 

Proposition 4.2 The time averaged consistency error of the Leray model and Leray- 
alpha model satisfy, for smooth u, 

< / |t (u,u)|g?x >< C8 < ||u|| 2 >^ min{< ||Vu|| 2 >^,5 < ||Au|| 2 >^}. 
Jn 

Proof By the Cauchy-Schwarz inequality 

/ |to(u, u)|dx < ||u — u||||u||. 
Jn 

Both filters satisfy 

||u-u|| < C<5min{||Vu||,<5||Au||}, 

e.g., for the Leray-alpha model C — 1. Integrate this in time, apply the temporal 
Cauchy-Schwarz inequality, take limits superior and the result follows since the 
generalized limit is below the limit superior. 

It is useful to estimate the RHS of these bounds in terms of the Reynolds number 
and (after time averaging) for a general weak solution of the Navier-Stokes equations. 
To do so a selection of the reference velocity must be made. In this generality the 
natural choice is 

U :=< — / |u(t,x)| 2 dx >5 . 
L A Jn 

With this choice we have 

If, 1 1 

< — / To(u, u) dx >< C5 < — 3-||Vu|| — r u > , or 
L A J Q L2 Li 

1 f I / M , 5 V .._ .,, 1 1 . . . ,n 1 

<js J ko(u,u)|dx >< C— < -^||Vu|| 2 >2< j^\\u\\ 2 >2 . 

Rewriting these in terms of the non-dimensionalized quantities gives 

1 f S i i 

< — / |r (u,u)|dx >< —= e^U^. 

Now, in turbulent flow typically e m ^C- (e.g., by dimensional analysis, Pope [POO] . 
and experiment, |S84j . [S98] . The estimate e < C^- has also ben proven directly 
from the Navier-Stokes equations, e.g., jCKGOl] . |CD92] . [DF02] . |Wang97| ). Thus 

we have 

< ij j[ Mu,u)|rfx >< c± 2 ui = 4rJuI 



Since both LHS and RHS are quadratic in U the most incisive form of this estimate 
is 

1 fl l \ I A 5Reh 

< fwyj ro(u ' u)|rfx> "Z77F- 

Related and more detailed estimates can be obtained in the case of homogeneous, 
isotropic turbulence the techniques introduced in [LL 06bj. 

Remark 4.1 The above presentation has been based on considering the Leray reg- 
ularization's solution as an approximation of the true (unfiltered) solution of the 
Navier-Stokes equations. If the filter is at least formally invertible, Geurts and Holm 
jGHWfy. IGHQpy have shown that by a change of variables a related Leray-regularized 
problem can be constructed whose solution is an approximation to the average ve- 
locity u. This change of variables alters the definition of the consistency error as 
follows. If u = A~ x u then u = Au . Since w~u(as5^0 ) w — A~ l w ~ u 
and the equation w is a Leray-regularization approximating the filtered solution of 
the Navier-Stokes equations. Changing variables byw = AA~ l w = Aw in the model 
gives (after simplification): 

w ( + V • (wAw) — i/Af + Vq = f in Q. 

To calculate this model's consistency error, rewrite the space filtered Navier-Stokes 
equations as 

u t + V • (uAu) - z/Au + Vp - f = V ■ [uAu - uu] . 

In this way we define the tensor 

To(u, u) := uAu — uu 

as the consistency error of the Leray and Leray-alpha model considered as an ap- 
proximation to u . However, this change of variables is an isomorphism and thus 
does not alter the final result on the model's consistency error. 

4.2 Consistency error of the Leray-deconvolution model 

We now consider the consistency error of the Leray-deconvolution model and show 
that the (asymptotic as 5 — >• ) consistency error is 0(5 2N+2 ). To identify the 
consistency error tensor, the Navier-Stokes equations is rearranged 

u t + D N (u) ■ Vu - i/Au + Vp - f = V • [D N (u)u - uu] in Q x (0, T). 

The LHS is the Leray-deconvolution model and the RHS is the residual of the true 
solution of the Navier-Stokes equations in the model. Thus, the consistency error 
tensor is 

Tat(u, u) := Dn(u)u — uu for AT = 0, 1, • • • . 

As in the Leray model, adapting the analysis in [LL06aj . [LL06bJ to the present 
case, the model error is driven by the model's consistency error t^(u, u) rather 



than V ■ 7jv(u, u). Since tn = (Dn(u) — u)u the consistency error of the Leray- 
deconvolution model is dominated by the deconvolution error. As before, there 
are three general weak solution, solutions with the regularity typically 

observed in homogeneous, isotropic turbulence (the deconvolution error is estimated 
in [LL06b] and the same estimates hold here) and, to assess accuracy on the large 
scales, very smooth solutions. In this case the deconvolution error is bounded in 
[DE06], |LL06bj and induces a high order consistency error bound, given next. 

Proposition 4.3 The consistency error of the N th Leray- deconvolution model is 
0(5 2N+2 ) ; it satisfies 

J Mu,u)|dx <5 2N+2 \\A N+1 (-5 2 A + l)-^u\\\\u\\, 

<C^ +2 ||A* +1 u||||ii||. 
Proof By the Cauchy-Schwarz inequality and Lemma 2.3 

J |rjv(u, u)\dx < ||u — -Djy(u)||||u|| 

< 5 2JV+2 ||A JV+1 (-5 2 A + l)-^ +1 )u| 1 1 |u| | . 

In the case of homogeneous, isotropic turbulence, precise estimates of time aver- 
aged consistency errors can be given following |LL06b] . 

5 Conclusions 

The Leray-regularization approach shows theoretical promise and has appealing sim- 
plicity. The tests of the Leray regularization in Geurts and Holm [GH03] have also 
been positive and the initial tests in |LMNR06] the higher order models (13. ip have 
shown deamatic improvements over the Leray model ( 11. ip . Extensive and systematic 
computational testing is the natural next step. However, there are also substantial 
theoretical questions open as well. The first is to develop a similarity theory for the 
models (paralleling the similarity theory in [Mus96j for the Smagorinsky model, in 
[FHTOlj for the alpha model and in [LN06a] . [LN06b] for other deconvlution mod- 
els) to better understand how the regularization in the model truncates solution 
scales. The second natural question is to study other filters and deconvolution op- 
erators. At this point, we believe that extension to many (but not all) other filters 
should present only technical problems. Extension to other deconvolution operators 
is important. Due to the many approaches to solution of ill-posed problems, this 
extension will most likely be done on a case by case basis. In practical computations, 
most likely both 5 and ./V will be varying. Thus, a detailed understanding of limiting 
behavior in both variables is also an important open problem. 
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